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Abstract. We present a general review of the dynamics of topological solitons in 1 and 2 
dimensions and then discuss some recent work on the scattering of various solitonic objects 
(such as kinks and breathers etc) on potential obstructions. 
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1. Introduction 

Topological solitons arise in many areas of applied mathematics and in the mathematical 
description of some processes in physics - for a good review see [Ij. 
Q^ I They arise in the mathematical description of objects that are localised in space as the energy 

density of these objects is nonzero only in a finite region; ie it is significantly nonzero in a small 
region and goes to zero, exponentially or as an inverse power, as one moves away from this 
region. 

Their stability is guaranteed by topological considerations, normally associated with the 
I topology of 5^ —>■ maps. 

' The simplest examples of such maps (in 1+1 and 2+1) dimensions involve, respectively, the 

Sine-Gordon kinks and the solitons based on sigma models. 

As is well known the Lagrangian density of the (1+1) dimensional sine-Gordon model is given 

by 

5l ; C = d^(t)- d^"^ - A sin^ ((/.). (1) 

Solutions of the equations of motion which follow from ([T]) are well known. They involve kinks 
and antikinks, which are topological solitons; breathers, which can be thought of as bound states 
of kinks and antikinks, further bound states of kinks and breathers, and many other solutions, 
less interesting from our point of view. 

The Lagrangian for the (2+1) dimensional sigma model is given by: 

C = d^$- 9^0-^5 [(9^0 • d^'^f - {d^$- dj){d^'$- 5>)] - V{$), (2) 

where for V we can take any 'simple' function of (p^, which vanishes as (j)^ = 1. So, in this work 
we take 

V{$)=^^{l-4>l). (3) 
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Furthermore, we require that the vector lies on the unit sphere hence (j)- (f) = 1. These last 
two conditions make the static finite energy solutions of the Euler Lagrange equations which 
follow from ([2]), the solitons, topological and based on the topology of 5^ — > 5^ maps. The 
requirement that the static energy is finite forces the field 03 to go to +1 at spatial infinity thus 
compactifying and allowing us to consider (f) on the extended plane R^yjoo topologically 
equivalent to S"^. 

The three terms in ^ are, from left to right, the pure sigma model, the Skyrme and the 
potential term. The last two terms are needed to stabilize the solitons. As long as they are 
small they have no influence on the topology of the fields. 

Incidentally, other potential terms ^ have also been studied [2]. The results, in these cases, 
are only slightly different but their generic features are the same. Hence in this work we discuss 
our results obtained for (l3|) . In 3 spatial dimensions we have skyrmions and monopoles but they 
will not be discussed here; the interested reader can find many interesting details in [1]. 



2. Dynamics 

The dynamics of the Sine Gordon kinks is well known: the kinks refiect from each other and, in 
fact, not much can happen as the motion is in one dimension. 

The dynamics of the two dimensional sigma model solitons can be either, relativistic, {i.e. 
based on the Lagrangian above) or based on the Landau-Lifshitz equation which is given by 

where L stands for the spatial part of £, i.e. of ([2]). 

The dynamics of both cases is very different. In the relativisitic case we have the familiar 90^^ 
scattering. Thus, when two solitons are sent towards each other head on (i.e. at a zero impact 
parameter), the system evolves in such a way that after the scattering we have two outgoing 
solitons which are moving in the direction perpendicular to the motion of the original solitons. 

This has been explained in many ways; the most compelling one involves the 
indistinguishability of solitons As the system of two solitons is described by a function 
which is symmetric with respect of the interchange of the positions of these solitons, the phase 
space of their relative position, is really described by mod a reflection with respect to the line 
joining their positions. Hence, effectively, the space is where Z2 describes this reflection, and 
so it is a cone. A straight line motion in this space through the vertex of the cone corresponds 
to a 90'' motion when viewed in the 'opened-up' R^. 

In the nonrelativistic case the situation is different as the equations involve only first order 
time derivatives. Hence the motion takes place in a lower dimensional phase space. 

This has been analysed in some detail by Papanicolaou and Tamaras [3j who showed that 
when one introduces, for a system of two solitons, r = (xi, 2:2) - a 2 dimensional vector describing 
their relative position, the corresponding relative momentum p satisfies 

Pi ~ acijXj. (5) 

Hence the equation motion is of the form 

d xi dxj , . 

resulting in a motion along a circle. 

The dynamics of solitons in 3 spatial dimensions is even more complicated ^; however, many 
aspects of it can be related to the dynamics in 2 dimensions. 



All this discussion concerned solitons moving in a free space, i.e. in a space with no potential 
obstructions. In the next section we summarise our recent results on the scattering of individual 
solitons on a spatial obstruction - either in the shape of a potential bump or a hole. 

3. Potential Obstruction 

To introduce a potential obstruction we note that the energy of the soliton field, strictly speaking, 
is never zero, even though it vanishes exponentially as we move away from soliton's position. 
Hence the obstruction has to be introduced in such a way that it does not change the "tail" of 
the soliton i.e. it has to vanish when, in one dimension, ^ = or vr or, in two dimensions, when 

3.1. Sigma Models 

Although a choice of obstruction is highly nonunique, in our studies, we added an extra term to 
the Lagrangian and we chose it of the form a(l — cf^) {i.e. like ([3]) with a nonzero in some region 
of X and y). We chose this term to be independent of y so that it represented an obstruction, 
located in some finite region of x resembling a trough in the "hole" case or a dam in the "barrier" 
case. Then sending the soliton from a point far away from this obstruction we can study its 
scattering on this obstruction. 

We have performed many numerical simulations of such systems, varying both the sign and 
value of a and the velocity of the incoming solitons. In the barrier case we have found that the 
scattering is very elastic with the soliton behaving like a point particle, depending on its kinetic 
energy - either going over the barrier or being reflected by it. Hence the velocity of the outgoing 
soliton is always very close, in magnitude, to the original velocity. 

For the hole case the situation was different. Depending on the value of the velocity solitons 
were either tranmitted or trapped and sometimes were even reflected. In fig 1. we present plots 
of the solitons' positions (in rc as y = 0) as a function of time seen in some simulations {i.e. 
started with different values of initial velocities). 
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Figure 1. Trajectory of a soliton during the scattering for a well of width L = 10 and depth 
0.2, {fi reduced from 0.25 to 0.2). Here Os = 0.25 and a) v = 0.0106, h) v = 0.0102, c) v = 0.012, 
respectively 

In addition to the trapping, which shows that, this time, the soliton does not behave like a 
point particle, we have also seen reflections! We have studied this process in great detail and in 
[1] we argued that the interaction of the soliton with the hole proceeds through the excitement of 
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Figure 2. Outgoing velocity of the kink as a function of its initial velocity 

the vibrational modes of the soliton. This has led us to look at other models and, in particular, 
at the Sine Gordon model in one dimension, which does not have any genuine vibrational modes. 
We discuss the results in this model in the next section. 

3.2. Sine Gordon Model 

There are two easy ways of introducing the obstruction into the Sine-Gordon model: either by 
making A in ([1]) position dependent [5] or by altering the basic metric [6]. Here we discuss the 
results reported in [5]; the results of [6] are qualitatively very similar. We restrict our discussion 
to the case of the potential hole. 

The results of our studies have shown that in the sine Gordon model, as in the sigma model 
in two dimensions, the solitons can get trapped, be transmitted and bounce back. The process 
is inelastic and depends on the initial condition and on the size and the depth of the hole. If the 
initial condition of the soliton corresponds to an exact Sine-Gordon kink moving with a given 
velocity then there is a well defined critical velocity above which the kink get transmitted (with 
a certain loss of energy). Below this critical velocity the kink can be trapped or reflected. The 
ranges of velocities, at which the kink is reflected are very narrow. The value of the critical 
velocity decreases with the decrease of the depth of the well and as the hole becomes narrower 
the number of velocity windows corresponding to reflections gets larger. These windows tend to 
be very narrow indeed. In fig. 2 we present a typical plot of the outgoing velocity as a function 
of the incoming velocity for the case of a relatively narrow well (in this case A in the well is 
lowered from 1 down to 0.8) and the hole is relatively narrow - i.e. a soliton fits in it about 3 
times) . 

Thus we note that in the Sine-Gordon model the solitons, which have no vibrational 
modes, can still be reflected by a potential hole. Of course, although the solitons have no 
vibrational modes they have pseudovibrational ones which can get excited and radiate (similar 
conclusions, in a different context, were reached by Romanczukiewicz [7J). An example of such 
a pseudovibrational mode is the mode which describes the variation of the slope of the kink. 




The usual kink solution is given by 



(j){x) = 2 arctan(exp(0(x — xo)), (7) 

where xq is the kink position and 6 is its slope. For ([7]) to be a solution of the equation of motion, 
which follows from ([1]), we need to put 9 and A equal. However, if we put 6 different from A 
we excite the mode which corresponds to the variation of 9. In fact, looking at the simulations 
in detail we have noticed that when the kink enters the hole for which A is different from 9, it 
automatically tries to adjust its slope and so it excites this mode. Of course, as soon as this 
mode is excited the kink begin to radiate and it is the interaction of this radiation with the kink 
itself which is reponsible for the final outcome of the scattering process. 

4. Breathers 

The Sine-Gordon model, in addition to the kink, also possesses many other solutions. Amongst 
them are breathers which are given by 

(t>{x,t) = 2 arctan ( ^M^i) — j ^ ^^g^ 
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As their energy is 16\/l — uj"^ they are often thought of as bound states of a kink and an 
antikink with the binding increasing as ^ 1. Hence it is interesting to see what happens when 
a breather is sent towards a hole. 

We have performed several numerical experiments of such systems and we are now preparing 
a long paper with a description of our results [5] • So here we will mention only some of them. 
As expected our results have shown that the breathers do get trapped, pass with, in general, 
a different uj (note that increasing lo releases some energy), or split with either a kink or an 
antikink being ejected from the hole. 

In fig 3 we present a couple of pictures showing a breather just before a scattering on a hole, 
and some time later, in the case when the interaction with the hole leads to the splitting of the 
breather. 

5. Further solutions (wobbles etc). 

However, the kinks and breathers are not the only finite energy solutions of the sine-Gordon 
model. In fact, the model possesses also solutions which describe bound states of kinks and 
breathers. One of such solutions, the 'wobble' was recently discussed in detail by Kalberman 

m- 

In his paper Kalberman presents an analytic form of this solution, shows that it describes a 
static kink in interaction with a breather, and then discusses some of its properties. 

Recently we have looked at more general solutions of the sine-Gordon model describing 
kinks and breathers. Our work was based on the Hirota method [TU] of deriving such solutions. 
Using this method we obtained 
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where 

Ts = 1k{x -VKt} + rjK IK = cosh ax vk = tanh ax, (10) 

p^ ^oM^B-aK}- cos 9 ^ = 2ArcTan 1^1!^. (n) 

coshjcts — cxKj + cost' sinhjas — axj 
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Figure 3. Trajectory of a breather sent with velocity v = 0.06587375 towards a well of width 
L = 10 and depth 0.8 (ie A = 0.8); t = 600 and t = 800, respectively 
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This corresponds to the most general kink/breather field configuration in which the kink and 
the breather have arbitrary velocities (and so move with respect to each other). If we now take 
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Hence our expression becomes 
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where u; is a frequency varying from —1 to 1. This agrees with the expression given by Kalberman 
which describes a stationary field configuration in which the kink and the breather sit on 'top 
of each other' and are not in relative motion. 

As is clear from (I13p the field configuration depends on one parameter (the frequency of the 
breather) and so, as we showed in [H], we have studied the stability of this field configuration. 
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Figure 4. Energy as a function of time as seen in a simulation started with A = 1.15 

To do this we calculated (j) and its time derivative from p3|) and then used them as the initial 
condition for our simulations. The results of our simulations were in complete agreement with the 
analytical expression thus showing that the solution is stable with respect to small perturbations 
(due to the discretisations). 

Next we studied the stability of the wobble with respect to larger perturbations. In particular 
we changed the original slope of the kink: (i.e. in the expression p3|) we replaced exp{x) by 
exp{Xx) where A 7^ 1). As discussed in [llj we performed several simulations with A ranging 
from 1.05 to 1.3. 

Our initial conditions corresponded to incorrect initial field configurations which added some 
extra energy to the system. The system was then unstable and so it evolved towards a stable 
wobble emitting some radiation which was sent out towards the boundaries of the grid. For A 
close to one - the perturbations were small - hence the system returned to its initial configuration 
(with A = 1.0). For larger values of the perturbation the system was more perturbed and often 
not only kept on sending out its excess of energy but also, at regular intervals, altered its 
frequency of oscillation (increasing it) which allowed it to send out even more radiation. In fig 
4. we present the plots of the time dependence of the total energy as seen in the simulation in 
which A was set at 1.15. 

Next we performed a series of simulations in which the initial ('wobble') configuration was 
sent towards a potential hole. As for a single breather [8] we have found many possibilities. We 
have also found that the hole can separate the breather from the kink (in one simulation we even 
saw the kink being trapped in the hole while the breather bounced off the kink and returned 
to where the wobble originally came from). Our results on this are still preliminary as we have 
decided to postpone the full study of this problem to some future work. 

We have also looked at field configurations involving a kink with more breathers. In particular 
in [11] we gave an explicit expression for such a configuration corresponding to a kink and two 
breathers. We have also verified that this configuration is stable. 



6. The energy 

As we mentioned before the energies of the kink or the breather configurations have very simple 
forms. This has led us to enquire whether this is also the case for the more general field 
configurations mentioned above. In fact, we managed to show this: the energies of all field 
configurations that can be derived by the Hirota method are additive as they are determined 
entirely by the asymptotic values of the Tj functions that arise in the construction of these field 
configurations. The Tj functions arise when one sets (p = 2 log ^ and then solves the relevant 
equations (more details in p3])- Any interested reader can find more detail in [llj . 
In the cases we have discussed here the Hamiltonian density is given by 

H = {dt(ff + {d^(ff + [sin((/.)]2, (14) 

and so the energy becomes 

/oo 
dxH = 2d^ (In To + Inn) |^=-oo (15) 
-oo 

and so is determined entirely by the asymptotic values of Tj functions. 

In this way we have managed to show that the energies of solutions we have considered in 
this talk are: 



(i) For the 1-kink: 

(ii) For the breather 

fiii) For the wobble 
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-2'wobble — 8 + 16 =. (18) 



(iv) For a solution involving the kink and two breathers mentioned above (with their velocities 
set to zero) 



^^kink+2breathcrs = 8 + 16^ I - Ujj + 16^ I - UJ^, (19) 

where uji = sin^j, i = 1,2. 
7. Perturbed Field Configurations 

Given that the sine-Gordon model possesses many solutions which resemble perturbed kinks 
{i.e. which are given by kinks and breathers) we have also tried to see what happens when 
one perturbs a kink and lets it evolve in time. Again we have looked at various perturbations, 
paying particular attention to configurations which were generated by adding to the kink an 
extra perturbation of the form 

5^{t = 0) = — ^— , 5^{t = 0) = — (20) 
cosh(|Ux) at cosh(i/x) 

In our simulations we used various values of A, B, fj, and u. In all cases the perturbation made 
the kink move and generated many moving breather-like configurations. To see what the system 
would finally settled at we absorbed the energy at the boundaries of the grid. This had the effect 
of slowing down the kink and also of absorbing and/or altering some breather-like structures. 



The process was very slow and the results were somewhat inconclusive. What we can say at 
this stage is that a general field configuration gradually splits into moving kinks and breathers, 
and some radiation, which quickly moves out to the boundaries. However, the resultant field 
configuration is metastable; it still radiates, albeit very slowly, and gradually evolves towards a 
field configuration involving mainly a kink. Whether at the end of its evolution we end up with 
a kink or a kink with some breathers is hard to determine. 

8. Simple models 

Here we mention briefly our simple models which partially explain what we have seen in our 
numerical simulations. First we consider the sigma model in 2 dimensions. In this case the 
soliton possesses many vibrational modes [12] and so we have constructed a model which treats 
the soliton as a system of four masses (connected to each other by springs) [4]. The system 
of four masses is sent towards the potential hole and then as it falls into it the masses begin 
to oscillate. These oscillations then model the soliton vibrations seen in full simulations. The 
energy is transferred to these oscillations and if the energy of the centre of mass is too low the 
system is trapped in the hole. Sometimes, when the system reaches one of the edges of the hole 
it happens to be in a state that allows the energy of the oscillations to get transferred back to 
the system as a whole (the energy of its centre of mass) and the soliton can come out. Whether 
this happens or not depends on the flow of the energy between the variational modes and on the 
kinetic and potential energy of the soliton. Hence, as we showed in the model reproduces 
quite well the main features of the scattering pattern seen in the full simulations. 

In the Sine-Gordon case we have looked at the old results on the scattering of kinks on point 
inpurities [13] (showing a similar trapping/transmission/refiection pattern) and their recent 
explanation by Goodman and collaborators [H]. In [H] the authors explained the observed 
results by invoking an interaction of the kink with the oscillation of the vacuum (around the 
impurity) which was described by a standing wave whose amplitude was a further degree of 
freedom (in addition to the position of the kink). The model reproduced all the features of 
the results of the original simulations reported in [T3] and so the two models discussed by us 
in [3] are based on the adaptation of the ideas of Goodman et al to our case. In both models 
we introduced degrees of freedom describing various standing waves in the hole (in one model 
the waves were restricted to the edges of the hole and in the other they described the global 
standing waves in the hole). The waves in both models were chosen somewhat arbitrarily as the 
idea was not to reproduce the pattern in any detail but just to check whether the mechanism of 
Goodman can be applied to our case too. 

In fact both models worked surprisingly well. They reproduced the pattern quite well, 
although the critical value of velocity was a little too high. Given that these models involved 
only very few (3 or 4) degrees of freedom we were very encouraged by these results; they require 
further work to understand better which modes are important and which are less so. 

9. Final comments and Conclusions. 

We started this talk with a very brief review of the topological solitons and of their dynamics. 
Then we reviewed the results of our studies of the scattering of topological solitons on a potential 
obstruction, of both a barrier and a hole-type. We finished by reporting some results on breathers 
and wobbles (bound states of breathers and kinks). 

Our results have shown that when a soliton was sent towards the barrier its behaviour 
resembled that of a point particle. Thus at low energies the soliton was refiected by the barrier 
and at higher energy it was transmitted. The scattering process was very elastic. During the 
scattering the kinetic energy of the soliton was gradually converted into the energy needed to 
'climb the barrier'. If the soliton had enough energy to get to the 'top' of the barrier then it 
was transmitted, otherwise it slid back regaining its kinetic energy. 



In the hole case, the situation was different. This time, as the soliton entered the hole it 
gained extra energy. Some of this energy was converted into the kinetic energy of the soliton, 
some was radiated away. So when the soliton tried to 'get out' of the hole it had less kinetic 
energy than at its entry and, when this energy was too low, it remained trapped in the hole. 
During the scattering process, like in the case of a barrier, the soliton's size changed and so it 
started oscillating. Afterwards, even when the soliton left the hole, its size continued to oscillate. 
Hence some energy was transferred to the oscillations resulting in the emission of radiation i.e. 
in an inelastic process. 

We have also looked at the scattering of breathers on potential holes. As breathers can 
be thought of as bound states of kinks and antikinks they can split leaving a trapped kink or 
an antikink in the hole and allowing its partner to escape either forwards or backwards. In 
addition, as the energy of the breather depends on the frequency of its oscillation ('breathing') 
this frequency can change as well. All such phenomena were seen in our simulations and we 
hope to report on this more fully in the near future [8j. 

We have also looked at some field configurations, which are solutions of the equations of 
motion and which describe bound states of kinks and breathers. As the energy of each breather 
depends on its frequency (and vanishes in the limit of this frequency going to 1) the extra 
energy, due to these extra breathers, does not have to be very large. The solutions appear to 
be stable and this stability is guaranteed by the integrability of the model. We have tested this 
numerically and have found that small perturbations, due to the discretisations, do not alter 
this stability. To change it we need something more drastic - like the absorption or the space 
variation of the potential {i.e. the coefficient of the sin? term in the Lagrangian). But even then 
the effects are not very large - one sees splitting of breathers etc but no 'global annihilation'. 

Finally, all our results (on the scattering of topological solitons on obstructions) generalise 
to other models; such as i.e. the Landau-Lifschitz model with a position dependent potential 
or an external magnetic field, other models in (1+1) dimensions, such as a X4>'^ model or even 
models describing ferro- and anti-ferro- magnets |15j . 

Clearly a lot of work still has to be done in this area. 
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